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INTRODUCTION 

The purpose of this paper is to recommend an algorithm to be used on 
data from a dual-frequency fluorosensor (i.e. one using two wavelengths for 
excitation of chlorophyll-^ fluorescence) to compute total chlorophyll-^ 
concentration and to partition that chlorophyll between two color groups 
present in a mixed phytoplankton population. The recommended algorithm is 
based on laboratory and field-testing experience gained with the Airborne 
Lidar Oceanographic Probing Experiment (ALOPE) fluorosensor at NASA's Langley 
Research Center. 

As with the single-laser fluorosensor, an assumption must be made that 
the fluorescence efficiency of each color group remains constant over the 
area being "calibrated**, but the dual-frequency technique can account for a 
shift in the overall or net fluorescence efficiency that would result from 
a shift in the relative abundance of two populations having different 
efficiencies. Therefore, the two-frequency technique can provide more 
accurate total chlorophyll estimates even if there is no interest in parti- 
tioning into color groups. 

Partitioning of the total chlorophyll into color group components requires 
knowledge of the spectral characteristics of the fluorescence excitation 
spectra of the color groups. Techniques used with single-laser f luorosensors 
to calibrate the fluorescence signal using concurrent measurements of 
chlorophyll ^ (sea truth) cannot be extended to the dual frequency technique 
by simply changing the dimensionality of the equations. This is because there 
are no reliable techniques for providing sea truth values for the chlorophyll 
a. concentrations of the two color groups. That is, if total chlorophyll is 
comprised of chlorophyll from two color groups with concentrations and C2 

such that the total chlorophyll concentration is -h C2 » sea truth 

does not exist for and C2 . Conventional techniques for obtaining 

total chlorophyll extract the pigment from all cells indiscriminately. Micro- 
scopic phytoplankton identifications or cell counts provide information on 
the relative abundance of the various species, and these can be classified . 
(hopefully) into color groups, but the information on cell size distributions 
and chlorophyll per cell needed to translate this into component chlorophyll 
concentrations is virtually nonexistent. Clearly it is beyond the capability 
of conventional shipboard techniques. 

This paper will describe algorithms for computing -, C^ 9 and C2 

at progressively more quantitative levels depending on the amount of infor- 
mation available or assumed. The first or least quantitative level is that 
of real-time data that can show, without sea truth, the relative variation 



of 9 ^2 slIotl^ the flight track. At a higher level, an algo- 

rithm is presented for mapping total chlorophyll using two fluorescences 
which is more accurate than one using a single fluorescence when varying 
mixtures of phytoplankton color groups are present. This requires sea truth 
on total chlorophyll, but no additional assumptions. Finally, given sea truth 
on total chlorophyll and assumptions about the fluorescence excitation 
characteristics of the two color groups of phytoplankton present, an algorithm 
is presented for computing and C 2 as well as • In both algorithms 

that use sea truth, the criterion chosen as the basis for deriving model 
parameters is to select those parameters that minimize the total squared error 
in chlorophyll ^ measured at the sea truth stations. That is, if Cx,k > 
k = l,...,m , are the sea truth measurements of total chlorophyll at m 
sea truth stations, and C-p , k = l,...,m , are estimates of total 
chlorophyll based on the algorithm used, then the algorithm parameters are 
selected to minimize 


m 

Total Error - ^ 

k=l 


( 1 ) 


This is the conventional unweighted least-squares solution where "error" is 
defined in terms of total chlorophyll measurements. 


GENERAL LINEAR MODEL: THEORY 

Fluorometric techniques for measuring chlorophyll ^ concentrations in 
living phytoplankton cells (in vivo) are based on the assumption that if the 
cells in a fixed volume of water are excited by light energy, the induced 
fluorescent energy emitted by the chlorophyll ^ per unit of excitation 
energy will be proportional to the molecular density or volumetric concen- 
tration of chlorophyll _a. The fluorescence of chlorophyll _a molecules is in 
a narrow spectral range centered at 685 nm. If the excitation source is 
effectively monochromatic, such as that provided by a laser, and if the 
excitation light frequency is varied, a fluorescence excitation spectrum is 
generated. Peaks in the spectrum correspond to absorbance peaks of auxilliary 
pigments present in the cell. Therefore, because the four major color groups 
of phytoplankton are characterized by the presence or absence of the auxilliary 
pigments, fluorescence excitation spectra can be used as a means of classify- 
ing the color groups. It was this fact that was the basis for the development 
of the ALOPE fluorosensor which utilized four distinct excitation wavelengths 
that were selected to discriminate among the four major color groups of 
phytoplankton. In practice a fluorosensor flown on an aircraft at altitudes 
generally around 150 m (500 ft) fires light pulses from a laser into the 
water and senses the returned fluorescence in several spectral bands. To 
apply fluorometric techniques, the chlorophyll fluorescence at 685 nm must 
be normalized to correct for variations in the laser's penetration depth 
along the track. The best- known technique for accomplishing this is to 
divide the 685-nm fluorescence by the Raman scattering produced by the laser 
which, in theory, is proportional to the number of water molecules accessed 



by the laser or, equivalently, the penetration depth of the laser energy. 

This technique also corrects for any variations in the output energy of the 
laser. 

In a dual-excitation frequency technique corrections must also be made 
for spectral differences in the light penetration and for any differences in 
the excitation energies of the two sources. All of these corrections are 
accomplished if the Raman signal from each laser is used to normalize 
fluorescences produced by that laser. It will be assumed in the remainder 
of this paper that fluorescence refers to a normalized fluorescence where 
variations in penetration depth and excitation energy are accounted for. 

The linear model used for fluorescence from a single-wavelength excitation 
is 


F = b + aC^ (2) 

where the fluorescence F is a linear function of the chlorophyll-a concen- 
tration C-p . The term b represents a background fluorescence which is 
not related to chlorophyll ^ and the slope a is related to the fluorescence 
efficiency of the chlorophyll The parameters b and a are assumed to 
remain constant over a defined area. Concurrent measurements of Cr^ (sea 
truth) are regressed against corresponding fluorescences F , and the slope 
a and intercept b are estimated. 


The model for a dual-frequency fluorosensor is 


^il^l + ^12^2 


(3) 


where fluorescence resulting from laser i , Fi , is linearly related to 
the chlorophyll-^ concentrations of the two color groups, Cp and C 2 » 
and to a background fluorescence induced by that laser. Again, the 

parameters b^ , a^p and a ^2 assumed to be constants over the calibra- 

tion area, and the parameters a^^^ and a ^2 » sometimes called cross sections 
are related to the fluorescence efficiencies of the two color groups at 
excitation frequency i. 


Experience with the ALOPE fluorosensor has shown that the ratio of 
fluorescences R = ¥ 2 /^-^ can be used as an indicator of the relative abun- 
dance of the two color groups. In cultures of a single species grown in the 
laboratory, where background fluorescences are assumed to be zero (i.e., 
bi = b 2 = 0) , fluorescence ratios are equivalent to the cross-section ratio 
azj/aij for color group j . The two excitation frequencies 454 nm and 
539 nm have been used extensively in ALOPE field tests to differentiate the 
golden-brown and green phytoplankton color groups, the two color groups 
commonly found in coastal and shelf waters. In the absence of any background 
fluorescence, the fluorescence ratio F (539) /F (454) is approximately equal 
to 1.0 for a golden-brown population and 0.3 for a green population. 
Figure 1 is a plot of this fluorescence ratio along a 220-km long flight 
track that began in the Chesapeake Bay where a golden-brown diatom was the 



predominant species and ended in clear offshore shelf waters where green 
species formed a significant component of the population. Note a shift in 
the ratio from values near 1.0 toward lower ratios around 0.45. 

Before proceeding to a discussion of the algorithms, which assume that 

and a^j are constants, some discussion of this assumption is in order. 

The assumption that the fluorescence efficiencies of the color groups remain 
constant over the calibration area is necessary, just as it is in shipboard 
techniques that use flow-through f luorometers. Sea truth measurements of 
chlorophyll ^ serve essentially the same purpose as the chlorophyll extractions 
that are made periodically to calibrate the continuous fluorometer record. 

A much more serious problem may be the assumption that background 
fluorescence is constant. Dissolved organic matter is known to have a broad- 
band emission spectrum that can significantly contribute to fluorescence at 
685 nm, and, depending on the excitation wavelength used, the Raman signal 
can overlap the chlorophyll ^ fluorescence. Probably the best method for 
removing a variable background signal, i.e. isolating the fluorescence due 
to chlorophyll is to have sufficient spectral resolution in the emission 
spectrum above and below the 685 nm band. Then, as illustrated in figure 2, 
a varying background can be estimated and removed. In all the algorithms 
discussed here, b]^ and l >2 are assumed to be constants that can be 
estimated and removed by various techniques. 


ALGORITHMS 

Two situations will be considered in the algorithms that are presented. 
First, no assumptions will be made about the fluorescence excitation spectra 
of the two color groups. In this situation the algorithm can provide total 
chlorophyll, either its relative variation (without sea truth) or absolute 
variation (with sea truth). Second, it will be assumed that the ratios 
Rj = ^2j/^lj both color groups (j = 1 and j = 2) are known constants. 

In this situation both relative and absolute estimates of and C 2 can 

be derived. 


Situation 1. Estimation of Total Chlorophyll ^ 

Without Assumptions About Cross-Section Ratios 

Without sea truth data, plots of and/or F 2 versus distance provide 

information about the relative concentration of chlorophyll ^ along the 
flight track. Since a shift in the relative abundance of the two color 
groups with different fluorescence efficiencies can affect the assumed 
linearity of fluorescence with respect to total chlorophyll, a plot of 
versus distance can be used to delineate the portions of the flight track 
over which linearity can be assumed. That is, either F^ or F 2 will be a 
valid relative measure of total chlorophyll over a region, provided ^2^^1 
is fairly constant. 
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Note that to compute fluorescence ratios the data must consist of 
fluorescence pairs representing the same locations along the flight track. 
Since the lasers are fired sequentially rather than simultaneously some 
preliminary interpolation is required to estimate the fluorescence return 
of the unfired laser to pair with a measured return. 

Once sea truth values for total chlorophyll are available, it is 
recommended that a multiple linear regression equation 

S ■ »0 

be derived rather than a simple linear regression of Crj, on or F 2 

alone. The reason for this is illustrated in figure 3. The symbols repre- 
sent a hypothetical situation in which sea truth values of total chlorophyll 
concentration Crj, = + C 2 are plotted against a fluorescence computed 

from the model 

= 1.0 + 1.2C^ + 0.30^ (5) 

Line I represents the linear relationship that would exist if the population 
were exclusively from color group 1, and line II that for color group 2. 

Any mixture of groups 1 and 2 would result in a point lying between the 
two lines. 


The cluster of points in figure 3 in the range 7 j< C-j* ;< 10 was sampled 
from a patch that was predominantly color group 2, whereas the cluster near 
C-p = 3 was from a patch which was predominantly group 1. Since the 
efficiency of group 2 is lower than that of group 1, the increase in fluor- 
escence resulting from a tripling of the chlorophyll concentration was 
offset almost entirely by an overall reduction in fluorescence efficiency. 
Clearly a regression line drawn through these data would be a poor represen- 
tation of the true relationship, particularly outside the range of the 
measured chlorophyll, or if a high chlorophyll patch of color group 1 or a 
low chlorophyll region of color group 2 were encountered. 

A multiple regression of C.p on Fp and F 2 would prevent such errors 
At each sea truth station, the relationship between Crj. and F^]^, ^2]a 

should be modeled as 


S,k ^0 




^2^2k 


( 6 ) 


for k = l,...,m . The least-squares solution for the coefficients that 

minimizes errors in C , (see equation (1) ) , is given by 

i ,K. 

3 = (7) 

where B = (3q* Bi, 32^^ 5 ^ m x 1 column vector of sea truth 

chlorophylls Crp ^ . More explicitly, the formulas for the regression 
coefficients are’ 



- %<^2> 

- oo') 


6, - 


T 

Sy^d - P„S 


6„ - 


Ct - - 6/2 


( 8 ) 


where and SQrj, are the mean and standard deviation, respectively, of 

the sea truth chlorophylls; and are the mean and standard 

deviation, respectively, of the corresponding sample; pQ is the linear 

correlation coefficient of F^ and ¥2 at sea truth stations; and is 

the linear correlation coefficient of F^^ and C-p at the truth stations. 


Situation 2. Estimating C-p, C^, and C 2 
With Assumptions About Cross-Section Ratios 

It is now assumed that the ratio of cross-sections known 

constants Rp and R 2 for color groups 1 and 2, respectively. The equations 
governing the fluorescences are now 


F 

F 


1 

2 


bi + + a^2C2 


^2 + ^1^11^1 + ^^2^2^ 


(9) 


At this point some further assumptions must be made regarding the background 
terms bp and b 2 • One cannot solve for app, ^12 » ^1’ ^2 

total chlorophyll sea truth data alone. One solution is to find a minimum 
Fp and F 2 in the entire data set, assume that chlorophyll is zero at this 
location, and set 


b^ = min F^ 

(10) 

= min F^ 


Another solution is to assume a fixed ratio Rq = t> 2 /bp based on known 
spectral characteristics of, say, dissolved organic material. 

Assume first that bp and bo are estimated by other means (e.g., 
equation (10)) and then subtracted from the fluorescences. The resulting 
chlorophyll ^ fluorescences would then be given by 
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(11) 


^1 “ ^ 11^1 \ 2^2 


^2 ^ 2 ^ 12^2 


Solving these equations for Cj^ and C 2 gives 


Cl = - 
1 a 


11 


^2^1 ~ ^2 

R 2 -Ri 


S = a 


12 


^2 “ 


R2 -Ri 


and summing these gives 


T a. , a 


11 


12 


where 


and 


^2^ - ^2 
^1 = R„ - R. 


^2 - 

^2 = R_ - R. 


( 12 ) 


(13) 


(14) 


Although sea truth values are required to derive estimates of a^j^ . and a-j ^2 
note that U-|^ and U 2 are relative measures of and C 2 , respectively, 

accurate to within scale factors 1/^11 ^/^12 • , plots of 

and U 2 along a flight track can be computed to provide the relative 
variation of and C 2 • The sum of and U 2 » however, is not 

necessarily a good measure of total chlorophyll because the scale factors 
may differ significantly in magnitude. Like or F 2 9 the sum ^2 

is a good relative measure of total chlorophyll only in regions where 
F 2 /F^ is fairly constant , indicating a constant ratio. 


Given sea truth Crj. , k 
and U 2 


1 , . . . , 


U 


Ik 


^2^1k ~ ^2k 

R 2 - Ri 


Ti = ^2k ~ \^lk 

2k - R, 


and corresponding values of 


(15) 
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the least-squares solution for 1/^11 ^/^12 

= [u^u]"^ 


‘11 


L^12 


( 16 ) 


where U is the in x 2 matrix whose kth row is ^2k^ * 

the m X 1 column vector of sea truth chlorophyll measurements. 


More explicitly, the solution for a^^j^ and a ^2 

IV - (zuj^u^)^ 

ZU^^EU2^ - (ZU^U^)^ 

ZC^U^ZU^^ - IC^V^IV^V^ 


(17) 


where all sums are over k = l,...,m . With these parameters, then, the 
Uj^ and U 2 values can be converted to Cj^ = Ui/^H SiXid C 2 = U2/a22 • 

The final case considered assumes Rq = t> 2 /bj^ to be a known constant. 
Here the fluorescences are 


^1 “ ^1 ^11^1 ^12^2 


^2 ^0^1 \^ii*"i ^2^12^2 


Solving for , ^2 , and Crp gives 


(18) 


and 




11 


Co = - 
2 a 


^1 ^2 ~ ^0 
^11 .^2 “ ^IJ 

^ Po -^ i ' 


12 “12 


R 2 -Ri 


(19) 




11 


^12 " 


( 20 ) 


where 6 q is the sum of the constant terms in equation (19) . 


424 



A multiple linear regression of sea truth chlorophyll values on 
and U 2 similar to that described above in equations ( 6 ) to ( 8 ) would 
provide solutions for Bq , 6 ]^ = and 82 = ^/^ 12 - 

Explicitly, the solutions are 


a 


11 


S 






«T - "I'ni - ®2/“l2 

”1 ■ iOJ 


0'“11 


where = (R^ - / (^2 " ^i linear correlation 

coefficients for and C-p at sea truth stations. 


( 21 ) 
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Figure !•- Fluorescence ratios for the March 17, 1980 James River/shelf mission. 



If fluorescence emission spectrum is measured at wavelengths 
Aj = 685 nm, then an estimate of = background is: 

bj. = 

^j+1 ■ \i-l 


Figure 2.- Use of neighboring bands in emission spectrum to estimate background. 
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